Curreem et al. Cell & Bioscience 201 1, 1:16 
http://www.cellandbioscience.eom/content/1/1/16 



Cell & Bioscience 



RESEARCH Open Access 



General metabolism of Laribacter hongkongensis: 
a genome-wide analysis 

Shirly 0 Curreem 1 , Jade L Teng 1 , Herman Tse 1,2,3,4 , Kwok-Yung Yuen 1,2,3,4 , Susanna K Lau 1 ' 2,3,4 * and 
Patrick C Woo 1 ' 2 ' 3 ' 4 ' 



Abstract 

Background: Laribacter hongkongensis is associated with community-acquired gastroenteritis and traveler's 
diarrhea. In this study, we performed an in-depth annotation of the genes and pathways of the general 
metabolism of L hongkongensis and correlated them with its phenotypic characteristics. 

Results: The L hongkongensis genome possesses the pentose phosphate and gluconeogenesis pathways and 
tricarboxylic acid and glyoxylate cycles, but incomplete Embden-Meyerhof-Parnas and Entner-Doudoroff pathways, 
in agreement with its asaccharolytic phenotype. It contains enzymes for biosynthesis and p-oxidation of saturated 
fatty acids, biosynthesis of all 20 universal amino acids and selenocysteine, the latter not observed in Neisseria 
gonorrhoeae, Neisseria meningitidis and Chromobacterium violaceum. The genome contains a variety of 
dehydrogenases, enabling it to utilize different substrates as electron donors. It encodes three terminal cytochrome 
oxidases for respiration using oxygen as the electron acceptor under aerobic and microaerophilic conditions and 
four reductases for respiration with alternative electron acceptors under anaerobic conditions. The presence of 
complete tetrathionate reductase operon may confer survival advantage in mammalian host in association with 
diarrhea. The genome contains CDSs for incorporating sulfur and nitrogen by sulfate assimilation, ammonia 
assimilation and nitrate reduction. The existence of both glutamate dehydrogenase and glutamine synthetase/ 
glutamate synthase pathways suggests an importance of ammonia metabolism in the living environments that it 
may encounter. 

Conclusions: The L hongkongensis genome possesses a variety of genes and pathways for carbohydrate, amino 
acid and lipid metabolism, respiratory chain and sulfur and nitrogen metabolism. These allow the bacterium to 
utilize various substrates for energy production and survive in different environmental niches. 



Background 

In 2001, Laribacter hongkongensis, a novel genus and 
species that belongs to the Neisseriaceae family of p- 
subclass of the Proteobacteria, was discovered from the 
blood and empyema pus of a patient with underlying 
alcoholic cirrhosis [1]. Subsequently, it was observed 
that L. hongkongensis was associated with freshwater fish 
borne community- acquired gastroenteritis and traveler's 
diarrhea in human [2-6]. In addition to its capability of 
living under both aerobic and anaerobic conditions and 
in the intestines of human, a variety of freshwater fish 
and frogs, it can also survive and replicate as a free liv- 
ing bacterium in water obtained from drinking water 
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reservoirs [6-10]. Despite its capability of survival in 
diverse environmental conditions, it does not metabolize 
any sugar tested [1,3,4,11]. 

In this article, we present an overview of the general 
metabolism of L. hongkongensis based on the informa- 
tion obtained from its genome analysis. The metabolic 
pathways of L. hongkongensis were also compared to 
those of Neisseria gonorrhoeae, Neisseria meningitidis, 
Chromobacterium violaceum, Escherichia coli and Cam- 
pylobacter jejuni. N gonorrhoeae, N meningitidis and C. 
violaceum are the other three bacterial species in the 
Neisseriaceae family of p -Proteobacteria with complete 
genome sequences available [12-15]. N gonorrhoeae and 
N meningitidis are strict aerobes that have stringent 
growth requirements and humans are their only known 
reservoir and host [16]. Conversely, C. violaceum is 
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facultative anaerobic, highly versatile in its metabolism, 
and can be found abundantly in multiple ecosystems, 
including water and soil; and in tropical and subtropical 
regions [17]. E. coli is the prototype Gram-negative bac- 
terium with its metabolic pathways dissected in the 
greatest detail. C. jejuni is another Gram-negative, S- 
shaped, motile, asaccharolytic bacillus associated with 
gastroenteritis [18]. 

Results and discussion 

The general metabolism discussed in this context mainly 
focuses on coding sequences (CDSs) that were classified 
into Cluster of Orthologous Groups (COG) functional 
categories of group C (energy production and conver- 
sion), G (carbohydrate transport and metabolism), E 
(amino acid transport and metabolism), I (lipid trans- 
port and metabolism) and P (inorganic ion transport 
and metabolism). Overall, the number of CDSs that 
were classified into these COGs was 191, 92, 253, 85 
and 155 respectively in L. hongkongensis genome. These 
numbers are higher than those in N. gonorrhoeae and N. 
meningitidis genomes (113, 61, 146, 44 and 96 CDSs 
respectively and 120, 67, 156, 52 and 96 CDSs respec- 
tively) but lower than those in C. violaceum genome 
(208, 203, 423, 139 and 226 CDSs respectively). This 
large number of genes in C. violaceum genome is in line 
with its ability to survive in a wide range of environ- 
ments, whereas the relatively smaller number of genes 
in the genomes of N. gonorrhoeae and N. meningitidis 
reflects their fastidious growth requirements and limited 
host ranges. Phylogenetic relationship of L. hongkongen- 
sis and other bacteria that are included for comparative 
analysis in this study based on 16S rRNA gene shows 
that L. hongkongensis is most similar to C. violaceum 
within the Neisseriaceae family with complete genome 
sequences available (Figure 1). This taxonomic closeness 
is also reflected in its gene contents, in which L. hon- 
gkongensis shares the highest percentage of CDSs with 



C. violaceum (64.1%) compared with other bacteria 
(Table 1). 

Figure 2 illustrates the deduced central metabolism of 
L. hongkongensis from the genomic data. The recon- 
structed metabolic pathways reconciled the results of 
previous physiological and biochemical studies in light 
of the metabolic capacity of this strain and showed 
remarkable ability to adapt to diverse environment. 
Comparison with other Neisseriaceae bacteria, prototype 
E. coli and C. jejuni brings to light the unique capabil- 
ities of L. hongkongensis. Utilization of compounds as 
carbon and energy sources was more restricted than in 
C. violaceum, but more extended than N. meningitidis 
and N. gonorrhoeae. 

Carbohydrate metabolism 

In most bacteria, glucose is metabolized through Emb- 
den-Meyerhof-Parnas (EMP) pathway, a pathway that is 
also highly conserved in eukaryotes and archaea. Apart 
from EMP pathway, some prokaryotes adopt Entner- 
Doudoroff (ED) pathway for sugar metabolism by prim- 
ing of 6-carbon sugars and subsequent cleavage into 
two 3-carbon intermediates by aldolase enzymes [19]. 
While some bacteria, such as E. coli and Pseudomonas 
aeruginosa, have genes for both EMP and ED pathways, 
others contain only either one of them, or each pathway 
is dedicated for the metabolism of specific 6-carbon 
sugar [19]. In L. hongkongensis, both pathways were 
incomplete with regards to the predicted genes present 
in the genome (Figure 2). The absence of hexokinase, 6- 
phosphofructokinase and pyruvate kinase, which are key 
enzymes of EMP pathway, suggested that the bacterium 
cannot catabolize glucose through EMP pathway, in 
agreement with its asaccharolytic phenotype. Similar 
observations are found in other asaccharolytic bacteria 
such as C. jejuni, Bordetella pertussis, Bordetella para- 
pertussis and Bordetella bronchiseptica, in which the 
genes encoding glucokinase and 6-phosphofructokinase 
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Figure 1 Phylogenetic relationship of 16S rRNA among L hongkongensis HLHK9, N. meningitidis MC58, N. gonorrhoeae FA 1090, C. 
violaceum ATCC 12472, E. coli K12 MG1655 and C jejuni NCTC 11168. The tree was inferred from 16S rRNA data by the neighbor-joining 
method. Bootstrap values were calculated from 1,000 trees. The scale bar indicates the estimated number of substitutions per 50 bases. Names 
and accession numbers are given as cited in the GenBank database. 
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Table 1 Gene contents in L hongkongensis, N. meningitidis, N. gonorrhoeae, C. violaceum, E. coli and C. jejuni 



Total number of CDSs 

With homologues in L. 
hongkongensis 

RNA gene 

rRNA genes 
5S rRNA 
16S rRNA 
23S rRNA 

tRNA genes 



L hongkongensis C. violaceum ATCC N. gonorrhoeae 
HLHK9 12472 FA 1090 



N. meningitidis 
MC58 



E coli K12 
MG1655 



NCTC 11168 



3235 



7 
7 
7 

78 



4407 

64.1% 
(2075/3235) 



2002 
39.6% 

(1282/3235) 



4 
4 
4 
55 



2143 

40.0% 
(1295/3235) 



4 
4 
4 
59 



4321 
50.4% 
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24.6% 
(795/3235) 
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Figure 2 Central metabolic pathways of L hongkongensis HLHK9 deduced from its genome sequence. The intermediary metabolic steps 
of tricarboxylic acid (TCA) cycle, gluconeogensis, glyoxylate cycle, pentose phosphate pathway, fatty acid and amino acid metabolism are shown. 
The function predictions were based on sequence similarity to proteins with experimentally verified functions. Individual pathways are shaded as 
indicated: yellow for gluconeogenesis, green for TCA cycle, and blue for pentose phosphate pathway. Arrows indicate the flow of intermediates. 
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are also absent from their genomes [15,20]. For N. 
gonorrhoeae and N. meningitidis, although they both 
metabolize a variety of sugars, their genomes do not 
contain 6-phosphofructokinase gene. Instead, they use 
ED pathway for metabolizing the 6-carbon sugars 
[12,21]. In the genome of L. hongkongensis, the gene 
that encodes one (2-keto-3-deoxy-6-phosphogluconate 
aldolase), but not the other (6-phosphogluconate dehy- 
dratase), of the two key enzymes of ED pathway is pre- 
sent, suggesting an incomplete pathway exists for 
metabolizing 6-carbon sugars through this pathway. 
With the inability to utilize exogenous carbohydrates, it 
is thus expected that gluconeogenesis would be essential 
to generate 6-carbon intermediates required for other 
biosynthetic pathways in L. hongkongensis, and the pre- 
sence of complete gluconeogenic genes is coherent with 
this prediction (Figure 2). 

In line with the asaccharolytic nature of L. hongkon- 
gensis and C. jejuni, both genomes do not contain genes 
that encode a complete phosphoenolpyruvate-depen- 
dent, sugar transporting phosphotransferase system 
(PTS), one of the major carbohydrate transport systems 
in bacteria. The genome of L. hongkongensis contains 
genes that encode enzyme I, phosphocarrier protein 
HPr, HPr kinase/phosphorylase HPrK, enzyme IIA Ntr 
and enzyme IIA fru , whereas none of the genes that 
encode any of the components of PTS is present in the 
genome of C. jejuni [15]. Similar to L. hongkongensis, N. 
meningitidis and N. gonorrhoeae also contain cytoplas- 
mic PTS protein homologues including EI, NPr (a var- 
iant of HPr) and HPrK together with incomplete 
enzyme II homologues although they were both sacchar- 
olytic. The presence of other sugar permeases in their 
genome suggests that sugars are probably transported 
through other pathways. In the genomes of C. violaceum 
and E. coli, genes encoding complete PTS protein 
homologues (enzyme I and HPr) and multiple types of 
enzyme II permeases are present, compatible with their 
ability to use a large number of sugars [13,14]. In gen- 
eral, enzyme I and HPr are energy coupling proteins of 
the system while complete enzyme II complexes, PTS 
permeases, consist of IIA, IIB and IIC domains (as well 
as IID domain in the family of mannose). Instead of car- 
bohydrate transport, some enzymes may be involved in 
regulation [22]. The presence of enzyme IIA Ntr and 
enzyme IIA fru in L. hongkongensis suggested that they 
are probably needed for the regulation of other bio- 
chemical pathways rather than sugar transport. 

Genes encoding enzymes of both non-oxidative and 
oxidative branches of pentose phosphate pathway 
(PPP) and TCA cycle can be found in the genome of 
L. hongkongensis, C. violaceum and E. coli, whereas C. 
jejuni lacks the oxidative branch of PPP [11-15]. 
Although genes that encode enzymes for PPP and 



most of TCA cycle are present in the genomes of N. 
meningitidis and N. gonorrhoeae, the gene for malate 
dehydrogenase is absent from both genomes. The pre- 
sence of flavin adenine dinucleotide (FAD)-dependent 
malate:quinine oxidoreductase may substitute for the 
reaction catalyzed by malate dehydrogenase by oxidiz- 
ing malate to oxaloacetate, yet experimental analysis 
awaits to confirm such prediction [21]. When carbon 
sources are metabolized to acetyl-coenzyme A (acetyl- 
CoA) instead of pyruvate or phosphoenolpyruvate, 
TCA cycle intermediates cannot be replenished 
through anaplerotic reactions. A complete pathway for 
glyoxylate cycle found in L. hongkongensis provides a 
simple and efficient strategy for it to convert acetyl- 
CoA into anaplerotic and gluconeogenic compounds 
and enables it to grow with C-2 compounds as the 
sole carbon source. With glyoxylate pathway, gluco- 
neogenesis, non-oxidative PPP and anaplerotic reac- 
tions, L. hongkongensis will therefore be able to 
generate different intermediates required for biosyn- 
thetic metabolism. The presence of critical anaplerotic 
enzymes, including phosphoenolpyruvate synthase and 
phosphoenolpyruvate carboxylase, which convert pyru- 
vate and oxaloacetate respectively to phosphoenolpyru- 
vate, together with malate oxidoreductase in L. 
hongkongensis reflects its potential to grow on carbon 
sources such as malate (Table 2). Malate, together with 
succinate and fumarate, are C 4 -dicarboxylates. They 
are intermediates in TCA cycle that can be utilized by 
bacteria such as E. coli, B. subtilus and many Rhizo- 
bium spp. as nonfermentable carbon and/or energy 
sources under aerobic or anaerobic conditions [23]. A 
number of C 4 -dicarboxylates can be found in nature 
such as fermentation products (e.g. succinate) and 
complexing agents (e.g. oxalate and malate) [24]. The 
presence of several C 4 -dicarboxylates transporters rein- 
forces the prediction of using C 4 -dicarboxylates as car- 
bon sources in L. hongkongensis. This is in line with 
our experiments on minimal medium for L. hongkon- 
gensis which showed that L-malate can be used as sole 
carbon source [11]. 

Amino acid metabolism 

Similar to E. coli and C. jejuni, but not N. gonorrhoeae, 
N. meningitidis and C. violaceum, the genome of L. hon- 
gkongensis contains enzymes for biosynthesis of all 20 
universally found amino acids and also selenocysteine. 
Selenocysteine (Sec) is an amino acid present in many 
species in three domains of life [25,26]. It is encoded by 
opal codon (UGA), of which selenoprotein mRNA car- 
ries a selenocysteine insertion sequence element imme- 
diately downstream to the selenocysteine-encoding UGA 
codon [27]. Biosynthesis and incorporation of selenocys- 
teine requires four genes, including selA (encoding Sec 
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Table 2 CDSs related to anaplerotic reactions in L hongkongensis, N. meningitidis, N. gonorrhoeae, C. violaceum, E. coli 
and C. jejuni 



Enzyme 


Gene 


L hongkongensis 
HLHK9 


C. violaceum 
ATCC 12472 


N. gonorrhoeae 
FA 1090 


N. meningitidis 
MC58 


E. coli K12 
MG1655 


C. jejuni 
NCTC 11168 


Phosphoenolpyruvate 
carboxylase 


ppc 


+ 


+ 


+ 


+ 


+ 




Phosphoenolpyruvate 
synthase 


pps 


+ 


+ 


+ 


+ 


+ 




Phosphoenolpyruvate 
carboxykinase 


pck 










+ 


+ 


Malate oxidoreductase 
(EC 1.1.1.38) 


maeA 






+ 


+ 


+ 


+ 


Malate oxidoreductase 
(EC 1.1.1.40) 


moeB 


+ 


+ 






+ 





synthase), selB (encoding Sec-specific elongation factor), 
selC (encoding tRNA Sec ) and selD (encoding selenopho- 
sphate synthetase), whereas ybbB (encoding tRNA 2- 
selenouridine synthase) is needed for utilization. All the 
genes required for biosynthesis and incorporation of 
selenocysteine are present in the genomes of about 20% 
of bacteria with complete genome sequence available, 
with a majority in Proteobacteria and Firmicutes [28]. 
While L. hongkongensis is predicted to synthesize seleno- 
cysteine, such ability is not found in C. violaceum, N. 
gonorrhoeae and N. meningitidis which also belong to P- 
Proteobacteria. Although the reason for the selective 
presence of selenocysteine in some organisms remains 
unknown, it is usually present in the active sites of pro- 
teins with redox functions [27]. From the genomic data, 
it is predicted that only the most commonly found sele- 
noprotein, alpha subunit of formate dehydrogenase 
encoded by fdoG, is present in L. hongkongensis using a 
selenoprotein prediction software bSECISearch [28,29]. 

Similar to C. jejuni, the genome of L. hongkongensis 
contains a number of proteases, peptidases and trans- 
porters for degradation and transport of proteins or 
peptides into oligopeptides or amino acids. The amino 
acids can then be degraded intracellularly by various 
enzymes (Table 3). These amino acids can be used as 
carbon and nitrogen source for the bacterium. The pro- 
ducts of the reactions can enter central metabolic path- 
ways including TCA cycle and gluconeogenesis pathway, 
and can be used for anabolic and catabolic purposes 
(Figure 2 and Table 3). 

As predicted from the genomic sequence, several 
genes involved in the biosynthesis (amino-acid acetyl- 
transferase gene argA: LHK 02338, LHK_02366; acetyl- 
glutamate kinase gene argB: LHK_02337, LHK_02829) 
and catabolism (whole operon of arginine deiminase 
arcBIAICID: LHK_02727-LHK_02734) of arginine are 
duplicated in L. hongkongensis, suggesting the impor- 
tance of arginine metabolism in the bacterium. 
Previously, we have shown that the two isoenzymes of 



Af-acetyl-L-glutamate kinase (NAGK) encoded by dupli- 
cated copies of argB, NAGK-20 and NAGK-37, which 
catalyze the key reaction in the 8 -step arginine biosynth- 
esis, gave differential expression pattern in a compara- 
tive proteomic study of L. hongkongensis growing at 37° 
C (human body temperature) and 20°C (freshwater habi- 
tat temperature) [11]. With NAGK-20 showing a higher 
expression at 20°C and NAKG-37 showing a higher 
expression at 37°C, kinetic analysis revealed that NAGK- 
20 also had a lower optimal temperature for enzymatic 
activities and was inhibited by arginine whereas NAGK- 
37 had a higher enzymatic activity with a higher optimal 
temperature and was insensitive to arginine inhibition 
[11]. These observations suggest the two isoenzymes are 
involved in temperature adaptation. Further investiga- 
tion into the functions of duplicated genes in the meta- 
bolic pathway of arginine should yield fruitful insights 
into the lifestyle of L. hongkongensis. 

Lipid metabolism 

Fatty acids are synthesized via repeated cycles of con- 
densation, dehydration and reduction of carbon-carbon 
bonds. While the majority of bacterial membranes are 
composed of saturated fatty acids, presence of unsatu- 
rated fatty acids (UFAs) can increase the fluidity of the 
membrane [30,31]. Similar to the genomes of the other 
five bacteria, the genome of L. hongkongensis contains 
all the enzymes for biosynthesis of saturated fatty acids 
(Table 4). As for UFA biosynthesis, there are two 
known sets of enzymes that operate by different 
mechanisms: aerobic route and anaerobic route [32]. In 
the aerobic route, UFAs are formed by the oxidation of 
saturated fatty acids catalyzed by acyl-ACP oxidase 
(fatty acid desaturase). In the anaerobic route, trans 
double bond is introduced to the acyl chain of 
ft-hydroxy-decanoyl-ACP and isomerized to ds-3-dece- 
noyl-ACP by bifunctional 3-hydroxydecanoyl-ACP dehy- 
dratase/isomerase (FabA) (encoded by fabA), followed 
by 3-ketoacyl-ACP synthase I (FabB) (encoded by fabB) 
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Table 3 Amino acid catabolism of L. hongkongensis HLHK9 deduced from its genome sequence 



Amino 


Pathway/enzyme(s) involved 


Gene number 


Intermediates/products formed 


Pathways that intermediates 


acid 








enter 


Serine 


L-serine dehydratase 


LHK_02265 


Pyruvate, ammonia 


TCA cycle or gluconeogenesis 


Aspartate 


Argininosuccinate synthase 


LHK_021 72 


Fumarate, oxaloacetate, ammonia, 


~rr~ a , 1 

ICA cycle 








ATP 






Argininosuccinate lyase 


LHK_03122 








L-aspartate oxidase 


LHK 00001 








Aspartate aminotransferase 


LHK_01340 






Glutamine 


Glutamine synthetase 


LHK_01876 


a-ketoglutarate, ammonia 


TCA cycle 




Glutamate dehydrogenase 


LHK_01886 






• 

Glycine 


Glycine cleavage system P-protein 


Lni\_Uz/ zz 


A rv^ rvM~\ n i 2 C(~\ MAP)I-I 
/AlTirTlOr lid, LU2, IMAUrl 


.... 
Ammonia assimilation 




Glycine cleavage system H protein 


LHK_02723 








Glycine cleavage system T protein 


LHK_02724 






Alanine 


Alanine dehydrogenase 


LHK_02210 


Pyruvate 


TCA cycle or gluconeogenesis 




Alanine racemase 


LHK_00350 








D-amino acid dehydrogenase 


LHK_00934 






Glutamate 


Glutamate dehydrogenase 


LHK_01886 


a-ketoglutarate, ammonia 


TCA cycle 


Arginine 


Arginine deiminase pathway 




Ammonia, ATP 


Ammonia assimilation 




Arginine deiminase 


LHK_02729, 










LHK_02734 








Ornithine carbamoyltransferase 


LHK_02728, 










LHK_02733 








Carbamate kinase 


LHK_02727, 










LHK_02732 








Arginine decarboxylase pathway 




Putrescine 


TCA cycle 




Arginine decarboxylase 


LHK_01034 








Agmatinase 


LHK_01 140 






Proline 


Proline dehydrogenase 


LHK_01861 


a-ketoglutarate, ammonia 


TCA cycle 




1-pyrroline carboxylate 


LHK_01861 








dehydrogenase 









which catalyzes the elongation of ds-3-decenoyl-ACP to 
form UFAs. This pathway is well-studied in E. coli 
[32,33]. The presence of either mechanism is sufficient 
for the biosynthesis of UFA. Unlike C. violaceum, which 
has the genes for the aerobic UFA biosynthetic pathway, 
the genome of L. hongkongensis does not contain the 
gene that encodes desaturase. As for the anaerobic 



route, fabA-fabB is generally restricted to genomes of Gl- 
and y-Proteobacteria. Therefore, it is not surprising that 
the genomes of L. hongkongensis, N. gonorrhoeae, N. 
meningitidis and C. violaceum and that of C. jejuni (e- 
Proteobacteria) do not contain fabA-fabB. Alternative 
pathways for anaerobic UFA biosynthesis in the absence 
of fadA and fadB have been reported in other bacteria 



Table 4 Comparison of metabolic pathways for fatty acid metabolism deduced from the genomes of L hongkongensis, 


N. meningitidis, N. gonorrhoeae, C. violaceum, E. coli and C jejuni 




Pathway 


L hongkongensis 


C. violaceum ATCC N. gonorrhoeae N. meningitidis E. coli K12 


C. jejuni NCTC 




HLHK9 


12472 FA 1090 MC58 MG1655 


11168 


Fatty acid biosynthesis 








Saturated fatty acid 


+ 


+ + + + 


+ 


Unsaturated fatty 




+ - - + 




acid 








Cyclopropane fatty 
acid 


+ 


+ - - + 


+ 


Fatty acid catabolism 








Saturated fatty acid 


+ 


+ - - + 




Unsaturated fatty 




+ - - + 




acid 
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such as Streptococcus pneumoniae and Enterococcus fae- 
calis, in which a trans-2 to ds-3-decenoyl-ACP isomer- 
ase (FabM) can introduce a cis double bond into the 
growing acyl chain in the former and the presence of 
homologues of FabZ and FabF that can function as 
FabA and FabB in the latter [34,35]. Other options such 
as dual functions of FabF involved in saturated fatty 
acid biosynthesis as FabF and UFA biosynthesis as FabB 
has also been reported in Clostridium acetobutylicium 
and Lactococcus lactis [36,37]. While L. hongkongensis 
contains neither of the alternative fabA, fabB and fabM 
genes in its genome, whether its FabF has dual function 
is yet to be determined. This phenomenon has also 
been observed in the genomes of N. gonorrhoeae and N. 
meningitidis. It has been suggested that there may be 
uncharacterized enzymes and pathways for unsaturated 
fatty acid biosynthesis in these bacteria [38]. Analysis of 
membrane phospholipids in L. hongkongensis will reveal 
the types of unsaturated fatty acids that it possesses and 
help to delineate the possible biosynthetic pathway. 

In addition, the presence of two homologues of cyclo- 
propane fatty- acyl-phospholipid synthases (LHK_01324 
and LHK_03103) (CFA synthase) suggested the possibi- 
lity for L. hongkongensis to synthesize cyclopropane fatty 
acid. Cyclopropane fatty acids are found in the bacterial 
membrane and are believed to be involved in acid resis- 
tance in E. coli and in association with virulence and 
persistence of Mycobacterium tuberculosis in host 
[39,40]. 

For the catabolism of fatty acids, genes that encode 
complete set of enzymes for p -oxidation pathway of 
saturated fatty acids are present in the genome of L. 
hongkongensis (Table 4). This suggests that L. hongkon- 
gensis may utilize these fatty acids as an energy source. 
The product of p -oxidation of saturated fatty acids, 
acetyl-CoA, will enter TCA cycle. On the other hand, 
similar to N. gonorrhoeae, N. meningitidis and C. jejuni, 
the genome of L. hongkongensis does not contain any of 
the two genes that encode enzymes for catabolism of 
unsaturated fatty acids (Table 4). 

Respiratory chain 

The constituents of the respiratory chain of L. hongkon- 
gensis dissected under aerobic and anaerobic growth 
conditions revealed the metabolic potential to obtain 
energy from various sources. The primary function of 
respiratory chain of bacteria is to produce ATP through 
the transport of electrons from various electron donors, 
usually intermediates of metabolic pathways linked by 
quinones or ubiquinones, to various electron acceptors. 
The proton gradient generated between the cytoplasm 
and periplasm is then dissipated for ATP synthesis by 
FiF 0 ATP synthase, which is predicted to be present in 
the genomes of L. hongkongensis and the other 5 species 



compared in this context [32]. The genome of L. hon- 
gkongensis contains a plethora of dehydrogenases, 
including NADH dehydrogenase I, Rnf type electron 
transport complex, succinate dehydrogenase, formate 
dehydrogenase, proline dehydrogenase, electron-trans- 
ferring flavoprotein dehydrogenase, and D-amino acid 
dehydrogenase (Table 5). A variety of substrates as elec- 
tron donors, such as NADH, succinate, formate, proline, 
acyl-CoA and D-amino acids appears to be utilized. The 
Rnf type electron transport complex was found in the 
genomes of L. hongkongensis and E. coli, but not those 
of N. gonorrhoeae, N. meningitidis, C. violaceum and C. 
jejuni. The complex was first discovered in Rhodobacter 
capsulatus in 1993 [41]. Although its role is not fully 
established, it was proposed to be an alternative enzyme 
to NADH dehydrogenase for utilizing NADH as the 
electron donor and had been found in many bacteria 
[42-44]. 

As for the final step of the respiratory chain, L. hon- 
gkongensis genome encodes three terminal cytochrome 
oxidases, namely type aa 3 oxidase (a haem-copper oxi- 
dase), type cbb 3 oxidase (another haem-copper oxidase) 
and type bd oxidase (a quinol oxidase). These three 
cytochrome oxidases are responsible for carrying out 
respiration using oxygen as the electron acceptor under 
aerobic conditions (type aa 3 oxidase) and conditions 
with reduced oxygen tension (type cbb 3 and type bd oxi- 
dases) (Table 5). Among the three cytochrome oxidases, 
type c&Z? 3 oxidase is the most ancient one and is present 
in almost all Proteobacteria except anaerobic 5-Proteo- 
bacteria [45]. It is characterized by the high affinity for 
oxygen which is needed in microaerobic environment. 
In addition to cytochrome oxidases, L. hongkongensis 
genome also encodes a number of reductases [fumarate 
reductase, nitrate reductase, dimethylsulfoxide (DMSO) 
reductase and tetrathionate reductase], which allow the 
bacterium to carry out respiration with several alterna- 
tive electron acceptors to oxygen (fumarate, nitrate, 
DMSO and tetrathionate) under anaerobic conditions 
(Table 5). The presence of versatile reductases may give 
clues to the living habitats of the bacterium where alter- 
native electron acceptors can be found in that environ- 
ment. The presence of DMSO reductase and 
tetrathionate reductase is of particular interest. DMSO 
is abundant in aquatic environments which can be pro- 
duced from the photochemical oxidation of dimethyl 
sulfide (DMS) and also from eukaryotic microplankton 
[46,47]. It may reflect the environment where L. hon- 
gkongensis is found. As for tetrathionate reductase, the 
presence of a complete ttr gene cluster suggested that 
tetrathionate can be used as an electron acceptor during 
anaerobic reduction of tetrathionate to thiosulfate. Nota- 
bly, L. hongkongensis is the only bacterium that contains 
the ttr operon among the family Neisseriaceae members 
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Table 5 Comparison of components of the respiratory chains deduced from the genomes of L hongkongensis, N. 
meningitidis, N. gonorrhoeae, C. violaceum, E. coli and C. jejuni 



Pathway/enzymes 



L hongkongensis C. violaceum 
HLHK9 ATCC 12472 



N. gonorrhoeae N. meningitidis E. coli K12 
FA 1090 MC58 MG1655 



C jejuni NCTC 
11168 



Electron donors 

NADH dehydrogenase 3 I 

Succinate dehydrogenase + 

Rnf-type electron transport + 
complex 

Formate dehydrogenase + 

Proline dehydrogenase + 

Electron-transferring flavoprotein + 
dehydrogenase 

D-amino acid dehydrogenase + 

Electron acceptors for aerobic 
respiration 

Cytochrome oxidase {oo 3 Xype) + 

Cytochrome oxidase {cbb 3 type) + 

Cytochrome oxidase {bd type) + 

Electron acceptors for anaerobic 
respiration 

Periplasmic nitrate reductase + 
Nap type 

Nitrate reductase Nar type 

Fumarate reductase + 

DMSO reductase + 

Tetrathionate reductase + 

FtFqATP synthase + 



I and I 
+ 



nopA only c 



I and II 
+ 



I and II 
+ 



I and II 

+ 

+ 

+ 
+ 
+ 



a Two types of NADH dehydrogenases (NDH) can function in the respiratory complex I in bacteria. Type I NDH is a proton pumping NADH:ubiquinone 
oxidoreductase that transfer electron from NADH to ubiquinone for generation of electrochemical proton gradient [75]. Type II NDH is a monomeric enzyme in 
which NADH oxidation is not coupled to proton pumping [76]. 

b The gene cluster that encodes NADH dehydrogenase I in most bacteria are consisted of 14 different genes. C. jejuni contains 12 of the 14 genes in which nuoE 
and nuoF are absent [15]. 

c The basic gene components for periplasmic Nap type nitrate reductase to be regarded as functional are napABCD [68]. C. violaceum contains only napA. 



with complete genome sequences available (LHK_01476- 
LHK_01478). Unlike the genes for carbohydrate, lipid 
and amino acid metabolism and those of the respiratory 
chain which are phylogenetically most related to other 
members of the Neisseriaceae, these three genes are 
most similar to their homologues in Thiobacillus deni- 
trificans ATCC 25259 and Aromatoleum aromaticum 
EbNl which were both P-Proteobacteria (Figure 3, 4 
and 5) [48,49]. The ability to respire tetrathionate is the 
characteristics of certain genera of Enter obacteriaceae 
including Salmonella, Citrobacter, Yersinia and Proteus 
[50]. It was recently reported that tetrathionate respira- 
tion provided growth benefit to S. Typhimurium in the 
lumen of inflamed intestine over other commensals/ 
microbiota inhabited in the gut [51]. Previously, only 
the source of thiosulfate, but not tetrathionate, was 
known in mammalian hosts. Toxic hydrogen sulphide 
gas produced by colonic bacteria in the intestine is 
detoxified by colonic muscosa into thiosulfates, thus a 
large quantity is present in the intestinal lumen [52,53]. 



It is now established that thiosulfate can be converted 
to tetrathionate upon reaction with reactive oxygen spe- 
cies that are generated during inflammation [51]. The 
presence of complete ttr gene cluster in L. hongkongensis 
genome suggests that it may confer survival advantage 
in mammalian host in association with diarrhea. On the 
other hand, the source of tetrathionate in nature is not 
clear, but it seems likely to occur in bacterial commu- 
nities that include sulfate-reducing bacteria as it has 
been detected in humid soils that support growth of 
such bacteria [50]. 

Sulfur metabolism 

Sulfur is a constituent of amino acids cysteine and 
methionine, and also found in various coenzymes and 
other metabolites. The major form of inorganic sulfur in 
nature, sulfate, is usually incorporated into bacteria as 
organic compound through the pathway of sulfate 
assimilation. With the ability to synthesize L-cysteine 
and L-methionine, L. hongkongensis is expected to 
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1 oo i Salmonella enterica subsp enterica serovar Typhimurium str. LT2 (NP_460348) 



100 



51 
100 



93 



97 



93 



94 



100 



100 



100 



95 



100 



100 



100 



89 



100 



87 



100 



L Salmonella enterica subsp enterica serovar Typhi str. Ty2 (NP_805054) 
Enterobacter cloacae SCF1 (YP_003941 684) 

- Proteus mirabilis HI4320 (YP_002151412) 

- Yersinia enterocolitica subsp enterocolitica 8081 (YP_00 1005907) 

- Serratia proteamacuians 568 (YP_00 1477301) 

Pasieureila multocida subsp multocida str. Pm70 (NP_245658) 

Haemophilus somnus 1 29PT (YP_71 91 07) 

1001— Shewanella baltica OS185 (YP_00 1367923) 

~~ I— Shewanella putrefaciens CN 32 (YP_001182191) 

Aeromonas hydrophila subsp hydrophila ATCC 7966 (YP_857776) 

Edwardsiella tarda EIB202 (YP_003295701) 

Photobacterium profundum SS9 (YP_1 31676) 

- Vibrio vulnificus YJ01 6 (NP_93501 3) 

- Vibrio parahaemolyticus RIMD 2210633 (NP_798393) 
Acidovorax sp. JS42 (YP_986661) 

Comamonas testosteroni CNB 2 (YP_003278276) 

- Delftia acidovorans SPH 1 (YP_001 565262) 

■ Pseudomonas stutzeri A1 501 (YP_001 1 741 58) 

- Achromobacter xylosoxidans A8 (YP_003978476) 
r Bordetella parapertussis 12822 (NP_883349) 



98 
64 



100L Bordetella bronchiseptica RB50 (NP_887789) 

Magnetospirillum Magneticum AMB 1 (YP_422649) 

Laribacter hongkongensis HLHK9 (YP_002795474) 

Thiobacilius denitrificans ATCC 25259 (YP_315499) 

Aromatoleum aromaticum EbN1 (YP_1 58751) 

Anaeromyxobacter dehaiogenans 2CP 1 (YP_002492207) 



100 



• Hydrogenobaculum sp Y04AAS1 (YP_0021 21 582) 
— Woiineila succinogenes DSM 1740 (NP_907142) 
Vibrio harveyi ATCC BAA-1 116 (YP_001 445772) 



Figure 3 Phylogenetic analysis of ttrA encoded protein in tetrathionate reductase gene cluster in L hongkongensis HLHK9. 

Phylogenetic tree showing the relationship of ttrA encoded proteins of L hongkongensis HLHK9 to other species with complete ttr gene cluster 
constructed by neighbor-joining method. Bootstrap values were calculated from 1000 trees. The scale bar indicates the estimated number of 
substitutions per 10 amino acids. All names and accession numbers are given as cited in the GenBank database. 



undergo sulfate assimilation with reduction of inorganic 
sulfate to sulfide and enter L-cysteine and L-methionine 
biosynthetic pathways. There are two pathways for 
intracellular sulfate to be reduced to sulfite and further 
to sulfide in bacteria, which differ in the form of acti- 
vated sulfate as input (Figure 6A). The first pathway, 
commonly but not exclusively found in enteric bacteria 
such as E. coli and traditionally assumed to be used by 
the majority of bacteria, converts adenosine 5'-adenylyl- 
sulfate (APS) to S'-phosphoadenylylsulfate (PAPS) via 
adenylylsulfate kinase CysC, where PAPS is eventually 
reduced to sulfite by PAPS reductase [54]. On the other 
hand, the second pathway which was originally 



identified in plants and subsequently found in other 
bacterial taxa including P. aeruginosa, other members of 
a-, P- and y-Proteobacteria and M. tuberculosis, reduces 
APS to sulfite directly by APS reductase [55,56]. Bacter- 
ial PAPS reductase and APS reductase are highly homo- 
logous and the major difference is the presence of 
conserved two-cysteine motifs in APS reductase [55-57]. 
Phylogenetic analysis showed that APS reductase from 
plant and many bacteria such as Pseudomonas, Rhizo- 
bium (Sinorhizobium) and M. tuberculosis clustered 
together [55]. The enzyme is also assumed to be typical 
in P-Proteobacteria as phylogenetic analysis of the gene 
components suggests that homologues of APS reductase 
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10Q | Salmonella enterica subs p. enterica serovar Typhimurium str. LT2 (NP_460350) 
' Salmonella enterica subsp. enterica serovar Typhi str. Ty2 (NP_805056) 

Enterobacter cloacae SCF1 (YP_003941686) 

Proteus mirabiiis HI4320 (YP_002151414) 

Serratia proteamaculans 568 (YP_001 477299) 

Yersinia enterocolitica subsp enterocolitica 8081 (YP_001 005905) 

Pasteureila multocida subsp multocida str. Pm70 (NP_245660) 

Haemophilus somnus 1 29PT (YP_71 91 09) 



10 or Bordetelfa bronchiseptica RB50 (NP_887791) 



99 



100 



80 



-Bordetella parapertussis 12822 (NP_883351) 

- Achromobacter xylosoxidans A8 (YP_003978474) 

Delftia acidovorans SPH-1 (YP_001 565260) 



74 



■ Pseudomonas stutzeri A1 501 (YP_001 1 741 60) 

Comamonas testosteroni CNB-2 (YP_003278278) 



86 



87 



100 



Acidovorax sp. JS42 (YP_986659) 

Edwardsieila tarda EIB202 (YP_003295699) 

— Photobacterium profundum SS9 (YP_1 31 679) 

- Vibrio parahaemolyticus RIMD 2210633 (NP_798391) 
-Vibrio vulnificus YJ016 (NP_935015) 

Aeromonas hydrophila subsp hydrophila ATCC 7966 (YP_857778) 
Shewaneiia baitica OS185 (YP_00 1367921) 
Shewanella putrefaciens CN-32 (YP_001 1 821 93) 
Anaeromyxobacter dehalogenans 2CP-1 (YP_002492205) 



— T 

100 1 



62 



99 



100 



94 



Laribacter hongkongensis HLHK9 (YP_002795476) 

Magnetospirillum magneticum AMB-1 (YP_422652) 

Thiobacillus denitrificans ATCC 25259 (YP_315497) 

■ Aromatoieum aromaticum EbN1 (YP_1 58749) 

Vrtbrio /?an/ey/ ATCC BAA-1 1 1 6 (YP_001 445774) 

Hydrogenobaculum sp. Y04AAS1 (YP_002121580) 

Wolinella succinogenes DSM 1740 (NP_907140) 



0.05 



Figure 4 Phylogenetic analysis of ffrB encoded protein in tetrathionate reductase gene cluster in L hongkongensis HLHK9. 

Phylogenetic tree showing the relationship of ttrB encoded proteins of L hongkongensis HLHK9 to other species with complete frr gene cluster 
constructed by neighbor-joining method. Bootstrap values were calculated from 1000 trees. The scale bar indicates the estimated number of 
substitutions per 20 amino acids. All names and accession numbers are given as cited in the GenBank database. 



with the characteristic two-cysteine motifs are preva- 
lently found in this group of bacteria [55]. The charac- 
teristic two-cysteine motifs, CCXXRKXXPL and 
SXGCXXCT, are found in the C termini of all APS 
reductases, but not in PAPS reductases [56]. The two- 
cysteine residues are reported to be involved for the 
binding to the iron- sulfur cluster in plant, P. aeruginosa 
and M. tuberculosis, which contribute to the difference 
in substrate specificity [58,59]. However, as this pathway 
has not been well characterized among P-Proteobacteria, 
the role of APS reductase in sulfate assimilation is 
uncertain in these bacterial taxa. As PAPS reductase has 



absolute substrate requirement, the product of the APS 
kinase, PAPS, is essential for E. coli to assimilate sulfate 
[54]. While the genomes of C. violaceum and E. coli are 
predicted to contain homologues of sulfate adenylyl- 
transferase, APS kinase, PAPS reductase and sulfite 
reductase, which can encode for a complete sulfate 
assimilation pathway via PAPS, L. hongkongensis and N. 
meningitidis seem to utilize the other pathway. Sequence 
analysis revealed that L. hongkongensis and N. meningiti- 
dis contain homologues of APS reductase with the char- 
acteristic two-cysteine motifs, and the absence of APS 
kinase homologue in the genomes suggested that the 
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Aeromonas hydrophila subsp hydrophila ATCC 7966 (YP_857777) 
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58 
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95 



00L Shewanella putrefaciens CN-32 (YP_001 182192) 
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- Thiobacillus denitrificans ATCC 25259 (YP_315498) 

• Anaeromyxobacter dehalogenans 2CP-1 (YP_002492206) 
Hydrogenobaculum sp. Y04AAS1 (YP_0021 21 581 ) 



100 1- 



Wolinella succinogenes DSM 1740 (NP_907141) 
Vibrio harveyi ATCC BAA-1 1 1 6 (YP_001 445773) 



0.2 

Figure 5 Phylogenetic analysis of ttrC encoded protein in tetrathionate reductase gene cluster in L hongkongensis HLHK9 

Phylogenetic tree showing the relationship of ttrC encoded proteins of L hongkongensis HLHK9 to other species with complete ttr gene cluster 
constructed by neighbor-joining method. Bootstrap values were calculated from 1000 trees. The scale bar indicates the estimated number of 
substitutions per 5 amino acids. All names and accession numbers are given as cited in the GenBank database. 



reduction of sulfate to sulfite in these bacteria may not 
involve PAPS intermediates (Figure 6B and Table 6). 
Although the predicted PAPS reductase in C. violaceum 
also has the characteristic two-cysteine motifs of APS 
reductase (Figure 6B), it remains uncertain whether this 
protein catalyzes PAPS or APS. 

Nitrogen metabolism 

Nitrogen is an essential component in organisms that 
constitutes vital biomolecules such as protein and 
nucleic acids. Bacteria can utilize different kinds of 
nitrogen compounds, including inorganic compounds 
such as ammonia, urea and nitrate; and organic complex 



compounds such as amino acids and nucleosides, as 
sources of cellular nitrogen. There are two general path- 
ways found widely in bacteria for ammonia assimilation. 
After transport across the cytoplasmic membrane into 
the bacterial cell, ammonia can be incorporated into 
glutamate or glutamine by glutamate dehydrogenase 
(GDH) or glutamine synthetase/glutamate synthase (GS/ 
GOG AT) pathway [32]. The two pathways differ accord- 
ing to the level of nitrogen present in the environment. 
When nitrogen concentration is high, ammonia assimi- 
lation via GDH is preferred as no ATP is consumed 
when glutamate is produced from oxoglutarate and 
ammonium. However, when nitrogen concentration is 
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PAPS reductase 
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Lh 


AQFGINGFYD 


SVAARQACCA 


ARKLEPLGRA 


LAGR- -GAWV 


TGLRRTQSV- 


TRSALAVRED 
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Ec 


EARYGKLWEQ 


GVEGIEKYND 


INKVEPMNRA 


LKELNAQTWF 


AGLRREQSG- 


SRANLPVLAI 
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Se 


EARYGKLWEQ 


GVEGIEKYNE 


INKVEPMNRA 


LKELNAQTWF 
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SRAHLPVLAI 
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Vc 


EARYGKLWEQ 


GVEGIERYNQ 


INKVEPMRRA 


LDELNISTWF 


SGLRREQSQ- 


SRASLPILSV 


175 


Sp 


EARFGQLWEQ 


GLDGLEQYNR 


LNKVEPMQRA 


LAELEVGTWF 


AGLRRSQAS- 


TREALPILAI 


189 


Sm 


AQYGMNGFYE 


SVEARHACCG 


VRKLKPLARA 


LDGA- -SYWI 


TGLRRGQSG- 


NRATTPFAEA 


174 


At 


RNKGLFSFYE 


- -DGHQECCR 


IRKVRPLRRA 


LKGL- -RAWI 


TGQRKDQSPG TCSEIPVVQV 


234 


Be 


AEHGLNAFYE 


SVELRKSCCH 


IRKVEPLNRA 


LADV- -GAWV 


TGQRREQSV- 


TRAELHEFEQ 


158 


Nm 


ESKGRFAFYD 


SVEARRECCR 


IRKIEPLNRA 


IAGA- -DAWL 


TGQRREQSA- 


TRTELPFAEY 


182 


Cv 


NFHGINGFYR 


SVELRKSCCA 


VRKLEPLKRA 


LAGR- -SAWI 


TGLRRQQSP- 


TRQDLADREY 


149 



Lh DPVNG LVKFSPL ADWSEAEVWE YVHRFSVPYN ALHDRHYPSI GCAPCTRAVA 201 

Ec QRGVF KVLPI IDWDNRTIYQ YLQKHGLKYH PLWDEGYLSV GDTHTTRKWE 222 

Se QRGVF KVLPI IDWDNRTVYQ YLQKHGLKYH PLWDQGYLSV GDTHTTRKWE 222 

Vc QNGVF KFLPV IDWTNKEVHY YLKDNDLPYH PLWEQGYLSV GDTHTTQKWQ 225 

Sp HGTRF KLLPI IEWSNKDVHL YLTQFNLPYH PLWDQGYVSV GDTHSSKPLE 239 

Sm DVERG LIKINPL ADWGIETIQA HVAAEGIPVN PLHSRGYPSI GCEPCTRAIK 226 

At DPVFEGLDGG VGSLVKWNPV ANVEGNDVWN FLRTMDVPVN TLHAAGYVSI GCEPCTRAVL 294 

Be DEARG IAKYNPL ADWTEADVWA YLKAFDVPVN PLHARGYPSI GCEPCTCAIR 210 

Nm DAGRG IGKYNPI FDWSEHDVWA YILANNVPYN DLYRQGFPSI GCDPCTRPVK 234 

Cv DADNG LVKFSPL LEWTEAEVWQ YLKENQVPYN ALHDRHYPSI GCAPCTRAIS 201 



Figure 6 Sulfate assimilation in L hongkongensis HLHK9. (A) Sulfate assimilation pathways in bacteria. Two pathways exist in bacteria which 
differ in the form of sulfate being activated. In one pathway that was once assumed to be the only one, APS is converted to PAPS by APS 
kinase and subsequently reduced to sulfite by PAPS reductase (green arrow). In the other pathway, APS is directly reduced to sulfite by APS 
reductase (orange arrow). Homologues of genes found in L hongkongensis HLHK9 were specified with gene numbers in blue. APS, adenosine 5'- 
adenylylsulfate; PAPS, 3'-phosphoadenylylsulfate. (B) Multiple sequence alignment of predicted APS reductase in L. hongkongensis with other 
characterized or related bacterial APS and PAPS reductases. Characterized APS reductases include those from Burkholderio cepacia, Arabidopsis 
thaliana, Rhizobium meliloti whereas characterized PAPS reductases include the one from £ coli. Only the selected region which contains two- 
cysteine motifs, CCXXRKXXPL and SXGCXXCT, in the C-terminal of APS reductase is shown. The conserved cysteine residues are shaded in yellow. 
The sequences were aligned with ClustalW. The abbreviations used and accession numbers (shown in parentheses) are as follows: At, A thaliana 
(GenBank:NP_1 93930); Be, B. cenocepacia J2315 [GenBank:YP_002231786]; Cv, C. violaceum ATCC 12472 [GenBank:NP_903244]; Ec, E. coli K-12 
MG1655, [GenBank:NP_417242]; Lh, L hongkongensis HLHK9 [GenBank:YP_00279461 1]; Nm, N. meningitidis MC58, [GenBank:NP_274183]; Se, 
Salmonella enterica subsp. enterica serovar Newport str. SL317 [GenBank:ZP_02697826]; Sp, Shewenella putrefaciens 200 [GenBank:ZP_01 705880]; 
Sm, Rhizobium {Sinorhizobium meliloti 1021) [GenBank:NP_385050] and Vc, Vibrio cholera TM 11079-80 [GenBank:ZP_04409338]. 
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Table 6 Comparison of sulfur metabolism deduced from the genomes of L. hongkongensis, N. meningitidis, N. 
gonorrhoeae, C. violaceum, E. coli and C. jejuni 



Pathway/enzyme 



Gene(s) 



L hongkongensis C violaceum 
HLHK9 ATCC 12472 



N. gonorrhoeae N. meningitidis E. coli K1 2 
FA 1090 MC58 MG1655 



C jejuni NCTC 
11168 



Sulfate assimilation 

sulfate cysD, cysN 

adenylyltransferase 

APS kinase cysC 

PAPS reductase cysti a 

APS reductase cysti a 

Sulfite reductase cysl, cysJ 

Sulfate transport 

ATP binding cysA 
protein CysA 

Permease protein cysil 
CysU 

Permease protein cysW 
CysW 

Sulfate binding cysP 
protein CysP 



a Bacterial PAPS reductase and APS reductase are highly homologous and both of their gene designation are cysti which may represent a PAPS reductase or APS 
reductase [58]. 

Predicted amino acid sequences encoded by cysH (CV_3574) in C. violaceum has the characteristic two-cysteine motifs of APS reductase and the presence of 
cysC suggests an ambiguity (represent by "?") in the identity of CV_3574 as PAPS or APS reductase. 



low, which is the case in most natural environment, 
ammonia is assimilated through GS/GOGAT pathway. 
L-glutamine is produced from L-glutamic acid and 
ammonia by GS in utilizing ATP and two molecules of 
L-glutamic acid are subsequently synthesized from L- 
glutamine and oxoglutarate by GOGAT because GDH 
has a lower affinity for ammonia. Similar to E. coli, L. 
hongkongensis has genes that participated in GDH and 
GS/GOGAT pathways, including gdhA that encodes 
NADP-dependent glutamate dehydrogenase, gltBlgltD 
that encode glutamate synthase and glnA that encodes 
glutamine synthetase (Table 7). On the other hand, the 
genomes of C. violaceum and C. jejuni contain only 
genes in GS/GOGAT pathway whereas those of N. 
meningitidis and N. gonorrhoeae contain solely genes of 
GS pathway. The ability of L. hongkongensis to utilize 
nitrogen is quite diverse. The existence of both path- 
ways in L. hongkongnesis suggests an importance of 
ammonia metabolism in the living environments that it 
may encounter. Ammonia can be found in natural 
environment such as natural water and being the form 
of excretory waste in fish, where L. hongkongensis had 
been reported to be isolated from [6,10,60]. 

In addition to its presence in the natural environment, 
ammonia can also be provided by alternative nitrogen 
sources such as urea. Urea can be hydrolyzed to ammo- 
nia and carbon dioxide by urease, in which a complete 
urease gene cluster is found in L. hongkongensis and 
none in N. meningitidis, N. gonorrhoeae, C. violaceum, 
E. coli and C. jejuni (Table 7). Ammonia formed in this 



process not only can provide nutrient nitrogen, but may 
also help L. hongkongensis to resist acid shock during its 
transit through the highly acidic environment of sto- 
mach by raising pH. With alternative way to metabolize 
different nitrogen source such as from urea, L. hongkon- 
gensis can gain some survival advantages under nitro- 
gen-limited conditions. The dual role in acid resistance 
and nitrogen metabolism of urease was evident in bac- 
teria such as Streptococcus salivarius, Helicboacter pylori 
and Yersinia enter ocolitica [61-63]. Investigation into 
the roles of urease in L. hongkongensis can provide 
insights towards the interplay of both functions and 
understanding about its survival strategy. Urea can be 
found in various environments in which L. hongkongen- 
sis can be found, including natural water and human 
host, where urea is present in a range of concentration 
in different parts ranging from saliva, stomach, blood to 
urine [64-67]. 

Alternatively, nitrate is another common form of 
nitrogen compounds found in the environment. By 
reducing nitrate, the majority of bacteria can incorporate 
nitrogen into building blocks, produce energy for cellu- 
lar processes or dissipate excess energy by respiration 
[68]. There are three types of bacterial nitrate reductases 
which are classified according to their localizations and 
functions, namely the assimilatory (Nas type), mem- 
brane-bound (Nar-type) and periplasmic (Nap-type), of 
which the Nas type is used for incorporating nitrate into 
building blocks [68,69]. In the assimilatory pathway, 
nitrate is first converted to nitrite by assimilatory nitrate 
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Table 7 Comparison of nitrogen metabolism deduced from the genomes of L hongkongensis, N. meningitidis, N. 
gonorrhoeae, C. violaceum, E. coli and C. jejuni 



Enzyme 


Gene(s) 


L hongkongensis 


C. violaceum 


N. gonorrhoeae 


N. meningitidis 


E. coli K12 


C. jejuni NCTC 






HLHK9 


ATCC 12472 


FA 1090 


MC58 


MG 1655 


11168 


Glutamate 


gdhA 


+ 




+ 


+ 


+ 




dehydrogenase 
















Glutamate synthase 


gltB/gltD 


+ 


+ 


- 


- 


+ 


+ 


Glutamine synthetase 


glnA 


+ 


+ 


+ 


+ 


+ 


+ 


Nitrate reductase 
















Membrane-bound 


narGHJI 9 




+ 






+ 




Nar type 
















Periplamsic Nap 


napABCD b 


+ 


nop A only 






+ 


+ 


type 
















Nitrite reductase 


nrfA/nifB/ 
aniA/nirK 




+ 


+ 


+ 


+ 




Nitric oxide reductase 


norB 




+ 


+ 


+ 






Urease 


ureABC 


+ 













a The basic gene components for membrane-bound type nitrate reductase to be regarded as functional [68] 
b The basic gene components for periplasmic type nitrate reductase to be regarded as functional [68] 
c The basic gene components for urease to be regarded as functional [77] 



reductase and then to ammonium by nitrite reductase 
such that it can be incorporated into cell materials. 
Absence of the genes that encode Nas type nitrate 
reductase and nitrite reductase in L. hongkongensis gen- 
ome suggested that no assimilatory pathway is present. 
For the respiratory pathway which mainly involved Nar 
and Nap type, while E. coli contains both types of 
nitrate reductases, L. hongkongensis and C. jejuni contain 
only periplasmic nitrate reductase whereas C. violaceum 
contains the membrane-bound type (Table 7). Even 
though C. violaceum contains napA that belongs to Nap 
system, the absence of other gene components suggested 
that this pathway is incomplete and probably non-func- 
tional. On the other hand, no homologues of nitrate 
reductases are found in N. meningitidis and N. gonor- 
rhoeae. The membrane-bound nitrate reductase, 
encoded by NarGHI operon, can be found in many 
nitrate-respiring and denitrifying bacteria and functions 
as the generator of proton motive force by coupling 
nitrate reduction in respiration. In contrast to Nar, the 
function of Nap is much diverse and appears to differ 
among bacteria. Apart from its involvement in anaerobic 
respiration in bacteria such as E. coli, it can have other 
physiological roles such as participating in redox balan- 
cing for optimal bacterial growth under certain physio- 
logical conditions and aerobic denitrification switch 
from aerobic respiration to denitrification or scavenges 
nitrate in some pathogenic bacteria [68]. The gene com- 
position and ordering of nap gene cluster shows hetero- 
genicity amongst different bacteria, with napD and 
napA as the only genes that are always found in the 
cluster [68]. In L. hongkongensis, it is predicted that nap 
operon contains 7 genes, napFDAGHBC, with the same 
gene composition and order as in E. coli. The ability for 



nitrate reduction in L. hongkongensis suggested that Nap 
type nitrate reductase is functional [1,3,11]. Apart from 
the above pathways, an absence of nitric oxide reductase 
in L. hongkongensis suggested that dissimilatory pathway 
for the production of dinitrogen is not operational. 

Conclusions 

The L. hongkongensis genome possesses a variety of 
genes and pathways for carbohydrate, amino acid and 
lipid metabolism, respiratory chain and sulfur and nitro- 
gen metabolism. These allow the bacterium to utilize 
various substrates for energy production and survive in 
different environmental niches. 

Methods 

CDSs identified in the L. hongkongensis genome were 
annotated as described in our previous publication and 
classified functionally according to the COG methodol- 
ogy [11]. Annotated genes were mapped to pathways 
according to the Kyoto Encyclopedia of Genes and Gen- 
omes (KEGG) database and MetaCyc to help identify 
metabolic pathways, and refined with thorough literature 
mining and experimental data [70,71]. When "gap" 
exists in the metabolic pathways, reciprocal-best-hit 
search was adopted to identify orthologues using experi- 
mentally verified protein sequence as query. CDSs 
belonging to metabolism-related COG clusters (C, G, E, 
I, and P) were selected for further examination and 
review. Other CDSs of potential interest were identified 
by keyword search using the names of common amino 
acids, carbohydrates, nucleotides, coenzymes, lipids, 
inorganic ions and metabolites. Manual annotation and 
analysis of the assigned function was performed by 
sequence similarity search using BLAST against the 
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NCBI nr database, and assisted by conserved domain 
search (CD-search), identification of signature sequence 
motifs and sequence analysis using InterProScan. Locali- 
zation patterns were predicted using PSORTb where 
appropriate [72]. Comparison of CDSs belonging to 
metabolism-related COG clusters among L. hongkongen- 
sis HLHK9, C. violaceum ATCC 12472, N. meningitidis 
MC58, N. gonorrhoeae FA 1090, E. coli K12 MG1655 
and C. jejuni NCTC 11168 was performed using KEGG 
database, MetaCyc and Integrated Microbial Genomes 
(IMG) database [73]. Phylogenetic trees were con- 
structed by the neighbor-joining method with MEGA 
4.0 using Kimura's two-parameter correction for 16S 
rRNA gene and Poisson correction for tetrathionate 
reductase [74]. 1553 nucleotide positions of 16S rRNA 
gene, and 1026, 233 and 389 amino acid positions of 
ttrA, ttrB and ttrC genes encoded proteins respectively 
were included in the analysis. 
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